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The dynamics of a two-mode Bose-Einstein condensate trapped in a double-well potential results 
approximately in an effective Rabi oscillation regime of exchange of population between both wells 
for sufficiently strong overlap between the modes functions. Facing this system as a temporal atomic 
beam splitter we show that this regime is optimal for a nondestructive atom-number measurement 
allowing an atomic homodyne detection, thus yielding indirect relative phase information about one 
of the two-mode condensates. 
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I. INTRODUCTION 

Since the first experimental achievements of Bose- 
Einstein condensation with a diluted atomic gas 
the investigation related to the detection of condensate 
relative phases and more specifically to the actual con- 
densate quantum state determination has generated a 
considerable literature 0, 0, 0] . A significant amount of 
discussion has been directed towards the detection of the 
relative phase of two Bose-Einstein Condensates (BECs), 
either in the form of interference between two indepen- 
dent BECs H 0, II H or due to indirect light probe of 
independent condensates T^. Essentially, acquiring in- 
formation about the BEC relative phase enables one to 
observe many interesting dynamical phenomena such as 
Josephson effect ^lIlllllfllMIIllIllIllIllllIM 

tand the transition from superfluid to Mott insulator 
. However the knowledge of such a phase could also 
be employed for partial or total inference on the conden- 
sate quantum state, through tomographic reconstruction 
as proposed by many authors in the past |3i 113 • 

Indeed, an atomic interferometric device implemented 
with high controllableparameters can be envisaged as an 
atomic beam-splitter 0,|2^|23|, which would be certainly 
important for schemes of quantum state reconstruction 
of BECs m m [23. It was previously noticed in Ref. 
[281 1 that a two-mode BEC trapped in a double-well po- 
tential could be envisaged as the atomic version of the 
Balanced Optical Homodyne Detection (BOHD), where 
the coherent tunneling of atoms would play the role of 
a temporal atomic beam-splitter. As it is well known, 
BOHD consists of mixing the signal field with a coherent 
local oscillator (LO) on a 50:50 beam-splitter to yield 
the necessary phase-sensitivity for signal field quadra- 
tures detection 30]. In the same sense, the phase 
of a signal atomic BEC mode could be thus determined 
by counting the difference of atoms in the wells of the 
trapping potential '2^. However in Ref. (2^ only an 
approximate description of a two-mode BEC was given, 
by neglecting cross-collision between atoms in different 
wells. As such, the approximate calculations were valid 
only for a small number of atoms (small condensate) and 
short time, determined by the ratio between tunneling 
and collision frequencies. Recently it was demonstrated 



that cross-collisions between atoms trapped in a double 
well potential could significantly increase the atom tun- 
neling rate for special trap configurations leading to an 
effective linear Rabi reg ime of population oscillation be- 
tween the trap wells [2lLl31| . This regime of oscillation is 
optimal for atomic homodyne detection of a signal BEC 
phase if the number of atoms in each well can be inferred 
from available experimental techniques. The Josephson 
coupling between distinct modes of a BEC trapped in a 
double well potential can be thus regarded as a reliable 
system for realizing a temporal atomic beam-splitter 
and consequently atomic homodyne detection. 

In this paper we develop the procedures for detecting 
a condensate relative phase through atomic homodyne 
detection. For the proper atomic homodyne detection 
process it is needed a secondary detection process able 
to determine the difference of atoms in the trap wells. 
The approach which we believe is the most promising 
in measuring the phase of a condensate is based on an 
extension of the homodyne measurements on a BEC pro- 
posed by Corney and Milburn in Ref. 0]. One of the 
wells of the double-well system is placed inside an opti- 
cal cavity, which is far off resonance with respect to any 
dipole transition in the atomic sample, allowing a dis- 
persive interaction between the light field and the atomic 
gas. Hence, the effect of the atoms is to shift the phase of 
the cavity field by a given amount dependent on the bal- 
ance of bosons in both wells, which may be measured by 
homodyne interferometry reflecting the internal dynam- 
ics of the condensate. Consequently, we simulate the ho- 
modyne current and its relation to the quadrature phase, 
possibly foreseeing experimental measured quantities and 
showing how the presence of cross-collisions enables a dy- 
namical regime ideal for such homodyne interferometry 
scheme. This paper is organized as follows. In Sec. II we 
review the main aspects presented in Ref. ^| for deriv- 
ing the effective stable Rabi regime. In Sec. HI we derive 
a detection model based on optical homodyne detection. 
In Sec. IV approximate solutions for the equations dis- 
cussed in Sec. HI are given, allowing the determination of 
the condensate relative phase through BOHD. In SEC. 
V the back-action on the condensate phase due to the 
continuous measurement process is analyzed and finally 
in Sec. VI a conclusion encloses the paper. 
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II. EFFECTIVE RABI REGIME AND BEG 
QUADRATURES DETECTION 

Before proceeding, we present briefly the quantum dy- 
namics of a BEC trapped in a double well potential as 
derived in [sij in order to justify the homodyne mea- 
surements discussed in later sections. The condensate 
model used in the follo wing discussion has been studied 
in previous papers llTl l3lj and so we only present 
an overview of it. Consider a double well potential trap- 
ping a Bose-Einstein condensate. The potential barrier 
is considered to be symmetric and the chemical potential 
is such that only two single-particle states are bellow the 
barrier separating the two wells, but in such a way that 
cross collisions between bosons of both wells may not be 
negligible. Those hypotheses enable a treatment of the 
many-body problem within a two-mode approximation. 
The well known bosonic many body Hamiltonian in the 
interaction picture is 



H 



2m 



*(r) 



j d3r*^(r)#t(r)<ir(r)*(r) 



(1) 



where m is the atomic mass, C/q 



measures the 



strength of the two-body interaction, a is the s-wave scat- 
tering length, and ^ are the Heisenberg picture field 
operators. If we consider a dilute gas in order that only 
s-scattering interactions are not negligible we may define 
the tunneling (fi), self-collision (k), and cross-collision 
(ry, A) rates, respectively as 



K 

A 



- I d^rul{r)[V{r) - V(^^\r - ri)]u2{r), (2) 



Uo 
2h 



(^) J d^ru*UiU*Uj, 



(3) 
(4) 
(5) 



where T^*-^^ is the harmonic approximation of the trap- 
ping potential around each minimum and Ui is the i-th 
mode function such that di{t) = J d'^ru*(r)'0(r, i). It is 
then possible to write down a two-mode single-particle 
Hamiltonian as 



H = ?l(2A(iV- 1) + Sl)[atS-|-6ta] + ?iry[at5. 



(6) 



where we have used the bosonic field operators relation 
to isomorphically map the many-body problem into the 
single-particle one. If we now introduce the Schwinger 
angular momentum representation 



1 



{h%-a)a). 



(7) 



Jy = ^(&^«-«^&), (8) 

J. - \{a^h^h'^a), (9) 
the Hamiltonian in Eq. (6) then becomes 

R^hyVL^ 2A{N - 1)]J, + AhrjJl + 2h{Ki - 77) (10) 

where we have neglected terms proportional to N and N"^ 
since they correspond only to a shift in the energy scale. 
The Casimir invariant is 



N ,N 



2 ' 2 



1), 



(11) 



which is analogous to an angular momentum model with 
total eigenvalue given by j = N/2. We may now use the 
Heisenberg picture to write the equations of motion for 
the angular momentum operators as follows 



-h^ + 2A(N - \)\Jy - AhrjiJy, J,l 



(12) 



Jy = h[n + 2A{N - 

4 = 2n{K- l])[jy,J^]+, 



2a(K-3r/)[J„ J,] + (13) 
(14) 



where [•, •]+ are anticommutators. This system of differ- 
ential equations can be solved numerically, and show a 
number of interesting elfects, such as self-trapping as dis- 
cussed in Ref. ISlI or in Refs. [Hill El Dj] 

in absence 

of cross-collision terms (77 = A = 0). Eqs. (12-14) es- 
sentially show that when the cross collision between the 
localized modes is taken into account the mode volume 
is increased. For a given fixed number of particles that 
means that the atomic density at each well is decreased 
and so the self-collisions rate occurring in each mode, as 
given by K — 77. Also the tunneling rate is increased as 
a consequence of the mode volume increase jsj as given 
by n + 2A{N — 1), and so dependent not only on the 
cross-coUisional rate A but also on the number of atoms 
in the trap. 

It is easily seen from the system of equations (12-14) 
that the presence of cross-collision inhibits self-trapping 
in the limit where k — 77 <C + 2A(A^ — 1), especially 
when rj ^ K. In such a case, Eqs. (12-14) result in 



I « 0, 



(15) 
(16) 
(17) 



being thus Jz approximately constant of motion. Here 
we have defined fl' = n + 2A{N - 1) + 8kJ^(0) as the 
new tunneling frequency, which explicitly depends on the 
cross-collision rate, on N, and on the initial condition 
for Jz- This regime could in principle be attained as 
discussed in Ref. 31] for special trap configurations. The 
new set of equations is thus easily solved to give 



X (t) = X (0) cos n't + jy{o) sin n't. 



(18) 
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If we suppose that initially both wells are equally popu- 
lated Ja;(0) = and the last equation reduces to 

X{t) =sm{n't)Jy{0). (19) 

As we shall see the dynamical regime imposed by Eq. 
(19) is optimal for atomic homodyne detection. 

By hypothesis we suppose that one of the two modes 
of the condensate (let us say mode B) is prepared in a 
coherent state H iM 113 in such way that /? = |/3|e*^ 
(Jy) can be rewritten as 

(Jy) = ^l/3| {{a^)e'' ' {&)e-'') - -|/3|(^«-./2>. (20) 

where Xg_T^/2 is the quadrature operator of the mode A. 
It is directly seen that 

iXit)) = \p\sm{n't){Xg_^/2), (21) 

which is the well-known result for balanced homodyne 
detection ^29^ ,30^ times a coherent amplitude dependent 
through Q' on the geometry of the trap, the total num- 
ber of particles and the initial condition of . It is also 
interesting to write the normalized operator Si = Ji/N 
with i G {x,y,z). The result in Eq. (21) means that 
even for large number of atoms, the self-trapping is to- 
tally suppressed and coherent oscillation takes place. We 
note that in this regime the frequency of oscillation in- 
creases with the total number of bosons in the system as 
n' = n + 2A{N -1) + 8NkSz{0), in such a way that the 
correspondent period decreases. This Rabi regime allows 
the double well trap to be envisaged as a realization of a 
temporal atomic beam splitter. Hence, the sine function 
modulating the homodyne current is an analogue to the 
beam-splitter transmissivity factor. For an ideal 50:50 
beam splitter the optimal situation would be such that 

n't^ {2n+l)n/2, (22) 

where n g TV. It is preferable to write Eq.(21) as 

- ^(4) (23) 

since the experimentally measured quantity would be the 
population difference given by the right-hand side of the 
above equation. Remark that since the above deriva- 
tion is for matter field instead of the BOHD and thus 
the quadratures are indeed given by combinations of the 
center of mass position and momentum operators for the 
mode A relatively to the mode B center of mass. Thus the 
atomic homodyne detection would essentially correspond 
to BEC mode A center of mass position and momentum 
measurements. 



III. ATOMIC HOMODYNE DETECTION 

For the complete implementation of the homodyne 
atomic detection a measuring process sensitive to the 



difference of atoms in the two modes is needed. Here 
we propose one possible implementation by letting one 
of the condensate modes to interact with a far off reso- 
nance cavity light field. The cavity output field is rccom- 
bined at a 50:50 beam splitter in a second stage BOHD 
as depicted in figure 1. The scheme is very similar to that 
proposed in Ref . jljj • One of the wells of the double- well 
system is placed in one arm of an optical cavity. The cav- 
ity is driven by a coherent field at the cavity frequency. 
It is supposed a dispersive interaction between the light 
field and the atomic gas in such a way that the field is far 
off resonance with respect to any dipole transition of the 
atomic species. Hence, the effect of the atoms is to shift 
the phase of the cavity field by a determined amount de- 
pendent on the balance of bosons in both wells. If the 
atom number in the cavity oscillates, so will the phase 
shift. Any tunneling of the condensate will be manifested 
in a modulated phase shift of the optical field exiting the 
cavity. In order to detect the light phase shift, it is con- 
sidered a common BOHD scheme. The light leaving the 
cavity is thus recombined with the reference beam in a 
50:50 beam-splitter and allowed to fall on the photode- 
tectors. 




FIG. 1; Atomic homodyne detection scheme. One of the BEC 
modes interacts dispersively with the cavity field, after [T^. 
The phase shift suffered by the cavity field is dependent on 
the imbalance of atoms in the BEC modes. Thus a secondary 
BOHD on the cavity output field allows inference of the con- 
densate relative phase. 

Throughout the following calculations we assume a 
bad-cavity limit, where any related atomic spontaneous 
emission rate is much smaller than the cavity field re- 
laxation rate, being thus neglected. In that regime the 
cavity field is undepleted and if the cavity light field is 
assumed to be far-detuned from any atomic resonance, 
the interaction Hamiltonian is effectively given [17| by 

H = J d\i>\v)[H,rr.-h^igir)c^c]i'{r) (24) 

where c and c' are the cavity field operators, g{r) is the 
intensity mode function and n = ^lj^/AA, with Rabi fre- 
quency flu, optical detuning A and Hem describing the 
center of mass motion. We may then write the above 
Hamiltonian in a single-particle formalism introducing 
the condensate operators a and and averaging over 



the optical mode function, resulting in the following in- 
teraction Hamiltonian 



(25) 



where ^ is the interaction strength. Since the cavity field 
is undepleted the total number of photons inside the cav- 
ity is a constant of motion. However the cavity field 
phase evolves with time. The phase time evolution can 
then be found by considering the Hcisenberg equation for 
the photon annihilation operator if the undepleted cav- 
ity field is assumed to be in a coherent state. Thus it is 
direct that 



(26) 



showing the direct dependence of the cavity field phase 
with the condensate imbalance operator (Jx)- If we now 
suppose the BEC is being monitored in a balanced homo- 
dyne way as in Fig 1 then it is a well known result that 
for balanced homodyne detection schemes, the difl[erence 
between both fields arriving at the photodetectors is pro- 
portional to the phase in such a way that 



iJ.f) = -\d\{X^-^/2) 



(27) 



where {Jxf) stands for the photon counting difference 
at the photodetectors, \d\ is the eigenvalue of the ref- 
erence beam annihilation operator, and X is the cavity 
field quadrature operator. In this last equation the light 
field phase (p varies with time depending on the conden- 
sate dynamics following Eq. (26). Hence, the measured 
photon difference gives us indirect information about the 
internal structure of the condensate since it relates it- 
self directly to the relative phase of the condensate in 
both wells of the trapping potential. A schematic circuit 
involving both the atomic and the optical homodyne de- 
tection process is depicted in Fig. 2. It is clear that (/> is 
a phase shift conditioned on the number of atoms in the 
BEC mode inside the cavity. 

We can develop this proposal further by writing 



(<^</.-7r/2(i)) = ('^<?i(t)-7r/2) 



4{t) 



(c)e- 



(t) 



(28) 

since the field that goes through the other branch of the 

beam-splitter docs not have a priori time dependence, 
being 0(i) as given by Eq. (26). If we assume both light 
beams being detected in a coherent state, then 

iXf) = -\c\\d\sm[ct>it)]- (29) 
We now have an explicit relation between the experi- 
mentally measured quantity {{Jxf)) and the condensate 
imbalance operator. In order now to access the relative 
condensate phase wc need to obtain a relationship be- 
tween the quadrature phase operator and the imbalance 
one. In the following we find such relationships and pro- 
vide some real insight on possible experiments measuring 
such quantities. 



FIG. 2: Circuit representing the atomic homodyne detection 
process. Grey linos arc for the condensate modes and dark 
lines for the optical fields. Beam-splitter 1 is the double- 
well trapping potential, while 2 and 3 are optical 50:50 beam- 
splitters. Atoms in the mode inside the cavity interact dis- 
persively with its field yielding a phase </> conditioned to the 
number of atoms in the mode. The light phase is detected 
through BOHD, thus giving information about the relative 
condensate phase. 



IV. APPROXIMATE SOLUTIONS 

Assuming again the cavity to be driven by a strong 
coherent field and being strongly damped, the cavity field 
is undepleted and in the k— eta <C limit the Heisenberg 
equations for the BEC operators are givenby 



-Vl'jy - 4k[Jj,, J^] + , 



(30) 

(31) 
(32) 



In the following we suppose that the BEC in the cavity 
is strongly embedded in the photon field in such a way 
that 



e = 



< 1 



(33) 



where A''^^ = {c?c). In other words we say that the density 
of the BEC is extremely small compared to that of the 
photons in the cavity. It is an experimental fact that such 
an approximation is quite correct since BEC densities 
range from 10^^ to 10^^ atoms/cm^. 

In the following, wc write down a solution for the set 
of differential equations given by Eqs. (21, 22, 23) up 
to first order in e by expanding the Schwinger operators 



clS J i 



(k) 



The zeroth order solution follows 



directly by integration and it reads 



{Jf\t)) = _|/3|(xf ^/^(t - 7r/2a;)), (34) 



where = Jl'^ -I- ^'^N'j. The first order solution may be 
found by considering the homogeneous solution (zeroth 
order) and applying the variation of parameters 



[1 -hcos2(2a;t)] 



3t 1 

! cos(2a;f) sin(2a;f) 

2 4a; 
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<V.> 

|c| |d| 




234 50 1 23 
tune time 

FIG. 3: Time evolution of the homodyne current as shows Eq. (43). For both graphics ^=0.01,f2'=25 Hz, uj—30 Hz and 
N=10000 atoms. In (a) we suppose a very large initial momentum (quadrature phase) of (Jy"''(0))=1667 and in (b) a small 
initial momentum of {J^°^(0))=0.001. 



. 3a; 



sin^(w<:) 



(35) 



Thus, up to first order in e the full solution reads 

(4)^(jW)+e(jW). (36) 

This is a complicated function of time but as it was al- 
ready mentioned, it shows that by mcasTiring the imbal- 
ance of population in the wells we acquire information 

I 



about the relative phase between both BEC modes as it 
is expressed in the quadrature phase operator -^6)_n-/2(i)- 
We now may plug those results in Eq. (26) and inte- 
grate it to obtain in zeroth order in e 



and in first order 



CO 



2 (^^-1/2(0) - f ) , (37) 



1 



4 ^ 16^2 



sin^(2a;t) 



. 3U! 
' o 



t sin(a;f) cos{wt) 



(^i-l/2W)}-(38) 



r 



Hence, up to first order the full solution reads 

(j) (j)^^"> + e(l)^^\ (39) 

As a matter of fact, it is possible to acquire information 
about the atomic gas quadrature from the quadrature of 
the light field by an atomic temporal homodyne scheme 
followed by a balanced light field homodyne detection. 
Though this solution is quite accurate, for real exper- 
imental data it may suffice to write down the solution 
up to zeroth order in e. Consider for instance k/O = 
0.02,^ = l{)-^Hz,n= 10^Hz,Nf = 10^°, n' = 9Hz and 
N=1000. For such a situation e 10~^ corroborating a 
zeroth order approximation of the problem. Henceforth, 
the phase (j) reads 



<^(i)^</)W(i) = C 



Nt 
2 



(40) 



It is interesting however to write down the imbalance in 
terms of the light phase since that phase could possibly 
measured by experimentalists giving indirect information 
about the atomic gas internal structure. In this sense it 
is easily seen that 



Nt 



(41) 



It is possible now to write down the full zeroth order 
solution to the Schwinger operators mean values so to 
write down an equation for the light field phase depend- 
ing explicitly on time: 

^cos(a;i)(Jr(0)) + ^), (42) 

meaning that for sufficiently large initial momentum 
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{Jy^\o)) the harmonic behavior should dominate and for 
smaller values the regime should be linear in time. 

For typical data of the system described here the zeroth 
order terms in e dominate and only coherent oscillations 
are observed for the imbalance of the atomic population 
given by the mean value of the Schwinger operator ■ In 
this approximation, the relationship between both phases 
is linear as we choose specific intervals of time when the 
potential barrier in the condensate acts as an ideal 50:50 
temporal atomic beam-splitter. These results make the 
model here presented an ideal system for optimal detec- 
tion of the condensed phase via two homodyne detec- 
tions: one temporal atomic one and a second one on the 
cavity output light field as depicted in Figs. 1 and 2. 

It is now possible to write down an expression which 
shows us how the light counting difference at the detec- 
tors is directly related to the internal structure of the 
condensate. Combining Eqs. (41) and (42) it is readily 
seen that 



iXf) = \c\\d\sm 



■ cos 



Nt 



The above expression shows clearly the relationship be- 
tween the measured homodyne current and the conden- 
sate quadrature phase. In this sense, it is possible to mea- 
sure the homodyne current {Jxf) and then acquire the 
necessary phase info about the condensate via Eq. (42). 
It is easy then to see that the evolution of the homodyne 
current with the condensate initial momentum (quadra- 
ture phase) is simply sinusoidal and its time evolution is 
quite similar but with the envelope function being mod- 
ulated as shown in Fig. 3 for two different choices of the 
initial momentum. The condensate relative phase is thus 
more evident for larger initial momenta. Such a result 
is evidently optimal for homodyne measurements and re- 
flects the effective Rabi dynamics discussed previously. 



V. MEASUREMENT BACK-ACTION 

We have considered a measurement process which al- 
lows the inference of a condensate relative phase through 
optical phase detection. However as is well known the 
continuous detection process induces a back-action into 
the condensate, altering thus the condensate phase dur- 
ing the measurement process. We now analyze those ef- 
fects by assuming the more realistic situation in that the 
cavity is driven by a strong coherent field of strength o, 
and is strongly damped at the rate 7. The procedure 
follows closely that by Corney and Milburn !17|. Hence 
the unconditioned evolution of the system (light field-t- 
BEC) is governed by the following master equation (tak- 
ing n = 1), 

hot = -«[^/,Pto*] +«C[c^cJj;,/5tot] - i(5 - ^)[c^c,/3tot] 



(44) 

where the initial detuning b ~ was chosen in order to 
remove the TV linear dependent dispersion. 




FIG. 4: Numerical calculation for the unconditional evolution 
of the homodyne current when F/il' — 0.0065, ^|cop/ri' = 
0.04. Time is normalized in units of fl' . 

It is possible to eliminate adiabatically the optical field 
from the master equation as in Refs. [tR Is^. under the 
hypotheses that the driving and damping terms dominate 
over the coupling one. This process leads to the master 
equation in terms of the atomic variables alone 



(45) 



where the measurement strength is F = 16f^<;^/7^ and 
I C|coKJn;> I — ^ I (large damping). We may then ob- 
serve the ensemble-averaged effect of the measurement 
in the operator moment equations for 77 ~ k, up to first 
order in e: 



Oy) 

(i) 



n'{X)+^\co\^{Jz) 



<\co\'{Jy)-^{J. 



(46) 
(47) 

(48) 



-ie\c' 



c, Ptotl 



C^Cptot - 0c) 



Those equations are the same for the case k — (co- 
herent oscillation) considered in Ref. but now the 
oscillation regime is attained in the presence of self- and 
cross-collisions, for 77 — > k and n — rj <^ fl' . We may now 
solve numerically this set of equations and with the aid 
of Eq.(26) and (29) find numerically the dependence of 
the homodyne current with the condensate phase quadra- 
ture. In Fig. 4 we depict the preselected homodyne 
current as a function of time for T/fi' = 0.0001 and 
rj/il' — 0.04. It is observed in general that even with the 
light field typical damping of the preselected state, the 
phase information of the field (relative to the condensate 
quadrature) is still present and shows similar behavior 
to that in Fig. (3a). The chosen initial momentum for 
the condensate is as such that the oscillatory behavior 
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shown in Fig. 4 is due both to the hnear term as 
well as the term proportional to the initial momentum in 
the sine argument of Eq. (43). 




FIG. 5: Numerical calculation for the conditional evolution 
of the homodyne current when T /fl' = 0.0001, rj/Q' = 0.04. 
Time is normalized in units of Q' . 

We now consider the postselected dynamics of the cav- 
ity field plus condensate system. It is usual then to nu- 
merically simulate stochastic realizations of quantum tra- 
jectories as already pointed out by several authors ^3^. Issl 
136 . 37, 38J. The resultant stochastic process is a diffusive 
evolution rather than the jump processes, which occur 
in the direct detection of atoms or individual photons 
since we have a condensate system continuously moni- 
tored by the optical homodyne detection scheme. In the 
presence of cross-collisions the conditional master equa- 
tion (which corresponds to an average over many runs of 
the experiment and many homodyne current records) for 
the optical field is given by 



§ ) = ^^M^^ 
dt J field 



V7^H[c]p., (49) 



where dW(t) is the infinitesimal Wiener increment '33, 
Pc is the density matrix that is conditioned on a particu- 
lar realization of the homodyne current u p to time t and 
V, Ti are the Wiseman's superoperators |33|- It follows 
then the conditional stochastic Schorodinger equation 



d\^cit)) = dt[-iH - -Jl 



/(t)Jx]|*c(i)>, (50) 



where H is given by Eq.(lO) plus Eq.(25) and ^c{t) de- 
scribes the conditional state of the system. The measured 
photocurrent is 



I{t) - 2T{X), + VvA(t), 



(51) 



where the stochastic term A{t) has the correlations 
{A{t)) = and {A{t),A{t')) = 6{t ~ t'). 

From Eq. (50) we see that the presence of cross- 
collisions in the rj {k — rj ^ Q') limit introduces 
only an harmonic correction due to the r]J^ term in 



the Hamiltonian. Then again, the cross-collisions repro- 
duce the homodyne interference pattern expected when 
there are no collisions at all {k — 0). Numerical simula- 
tions then show similar results to that found by Corney 
and Milburn |l7j for only coherent oscillation dynamics 
(k = 0) corroborating the fact that such system in the 
effective Rabi regime is optimal for homodyne detection 
in the sense that it attains a purely coherent oscillation of 
population dynamics without collapse and revival. The 
numerical simulations depicted in Fig 5 shows as well, 
how the experimental measured quantity given by Eq. 
(29) evolves in time under the conditioned detection for 
r/n' = 0.0001 and rj/n' = 0.04. We then observe that 
the homodyne current changes considerably in relation to 
the oscillatory unconditioned evolution. This is expected 
since the measurement back action alters considerably 
the whole system state. An efficient detection process 
with a larger initial condensate momentum would allow 
the evidence of the condensate phase. It may then be 
possible to experimentally confirm such results by mea- 
suring the homodyne current and then inverting Eq. (29) 
in order to access the imbalance of population between 
both wells and then to acquire the desired information 
about the relative quadrature phase of the Bose-Einstein 
condensates. 



VI. CONCLUDING REMARKS AND GENERAL 
DISCUSSION 

We have shown that in the effective Rabi regime of 
a double-well atomic BEC j^J an optimal condition for 
atomic homodyne detection scheme is found, which gives 
indirect measurement of the condensate relative phase. 
The double-well potential barrier acts as a temporal 
atomic beam-splitter with the transmissivity factor vary- 
ing with time and depending directly on the total number 
of bosons and cross-collisions strength by the corrected 
frequency fi' [23|. Up to first order in e, the Heisenberg 
equations of motion for the mean values of the Schwinger 
operators are exactly soluble even when the interaction 
with the light field is considered strong for sufficiently 
strong light intensities. 

Typical experimental data show that it may suffice to 
consider only zeroth order terms in the calculations which 
result in a linear relationship between the light phase 
and the condensate quadrature. In this sense it is sup- 
posed that the light phase may be detected with the aid 
of the scheme proposed in Figs. 1 and 2. It consists 
of a two stage homodyne detection, one optical and the 
other on the state of one of the two-mode condensate. 
Hence, we believe that such a system in this dynamical 
regime (effective Rabi) might be the appropriate choice 
to indirectly detect the relative phase between the two 
modes of a BEC in a double well potential in the form 
of Josephson-like tunneling in a regime of purely coher- 
ent exchange of population between both wells due to 
the strong presence of cross-collisions. Such a conclusion 
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is strongly supported by the calculations (analytical and 
numerical) discussed in this paper. 

Recently an outstanding experiment was realized based 
on stimulated light scattering to continuously sample the 
relative phase of two spatially separated atomic BECs 
that never interact • Our proposal on the other hand 
imposes that the two atomic BEC modes must be over- 
lapping in order that the effective (stable) Rabi regime 
(k — ry ^ ri') be attained. In face of Eq. (26) and (29), we 
expect experimentalists to be able to measure the relative 
phase with present technologies on trapping potentials. 
Experiments should then be able to detect the relative 
condensate phase, as given by the model here presented, 
possibly opening new frontiers in quantum phase engi- 
neering. We expect that these results may be useful in 
further experimental and theoretical studies on the state 



of a Bose-Einstein condensate as well as to applications 
on atom optics when such systems are extrapolated to 
an array of BECs. As a last comment, our results are 
certainly relevant for reconstruction and measurement of 
atomic quantum states psl |3 IHl may be useful 
for future implementations on quantum communication 
protocols pLIsgji. 
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